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We establish the existence of a chiral spin liquid (CSL) as the exact ground state of the Kitaev 
model on a decorated honeycomb lattice, which is obtained by replacing each site in the familiar 
honeycomb lattice with a triangle. The CSL state spontaneously breaks time reversal symmetry 
but preserves other symmetries. There are two topologically distinct CSLs separated by a quantum 
critical point. Interestingly, vortex excitations in the topologically nontrivial (Chern number ±1) 
CSL obey non-Abehan statistics. 



There has been enormous interest recently in the no- 
tion of spin-liquid states with fractionalized quasiparti- 
cle excitations, especially those that exhibit fractional 
statistics. Since, to date, unambiguous experimental ev- 
idence of fractionalization has only been achieved in the 
fractional quantum Hall effect (FQHE), it is reasonable 
to seek a high standard of theoretical control in mak- 
ing statements about the behavior of model systems that 
are being studied in this context. A breakthrough oc- 
curred when Moessner and Sondhi [l[ demonstrated the 
existence of a symmetry preserving spin-liquid phase in 
the quantum dimer model analogous to the short- 
range version Q of the RVB state proposed [H to be the 
key to understanding high temperature superconductiv- 
ity. This state has fractionalized excitations, but only 
with Fermi or Bose statistics. Motivated by a formal 
analogy between frustrated interactions and a magnetic 
field, Kalmeyer and Laughlin (KL) [f| proposed a varia- 
tional spin-liquid state for the triangular lattice quantum 
antifcrromagnct (AF) constructed in analogy with the 
Laughlin state 0] of the FQHE. This state breaks time re- 
versal symmetry (TRS) Q, and is thus a first example of 
a proposed class of chiral spin liquids (CSL) However, 
both analytical and numerical studies [§] suggest that the 
AF Heisenberg model on a triangular lattice possesses 
\/3 x a/3 Neel order, which breaks translational and spin 
rotational symmetries in addition to TRS. Therefore, it 
remains conceptually interesting to determine whether a 
CSL could exist as a stable ground state phase of any 
Hamiltonian with only short-range interactions. 

Here, we study the Kitaev model on the "triangle- 
honeycomb lattice" (Fig. HJa)) and show that its exact 
ground state is a CSL - i.e. a state which spontaneously 
breaks time-reversal symmetry but neither spin rota- 
tional nor translational symmetry, and which has frac- 
tionalized excitations. Indeed, there are two topologi- 
cally distinct CSL phases (with even and odd Chern num- 
ber, respectively), separated by a quantum phase transi- 
tion. For the topologically non-trivial (trivial) CSL, the 
vortex excitations obey non-Abelian (Abelian) statistics. 
(There are, however, an even number of electrons per 
unit cell, so this state may not deserve to be called a 
"Mott insulator".) As far as we know, this is the first 



model whose exact ground state is a CSL. 

The "triangle-honeycomb" lattice shown in Fig. [IJa) 
is obtained by replacing each site of a honeycomb lattice 
with a triangle - there arc thus 6 sites per unit cell. It 
is important to note that, in contrast to the honeycomb 
lattice, the elementary plaquettcs of this lattice have an 
odd number of sites. The Kitaev model (Eq. {]])) on this 
lattice is time reversal invariant. It is also invariant un- 
der inversion of the lattice and under 180° rotation of 
the spins along the x, y or z-directions, i.e. spin-orbit 
coupling has reduced the spin symmetry from SU (2) to 
D 2 . We will show that the ground state of the Kitaev 
model on this lattice does not break the spin symme- 
tries, nor any lattice translational or point group sym- 
metries. (Indeed, the spin correlation function is exactly 



zero beyond nearest neighbors [ill].) However, in both 
the topologically trivial and non-trivial states, the exact 
ground state spontaneously breaks TRS. Similar states 
were studied previously by Kitaev [l(j and Lee et al [l2| 
on a honeycomb lattice, but in models in which TRS is 
explicitly broken, e.g. by an external magnetic field. 

The existence of an exactly solvable model with a CSL 
ground-state is suggestive that such phases might exist 
in nature. One might try to realize something close to 
this precise model in a system of cold atoms in an optical 
lattice More generally, the fact that quantum Hall 
physics and non-Abelian anyons can occur without the 
need for a large applied magnetic field, opens the pos- 




FIG. 1: (a) The triangle-honeycomb lattice is constructed 
from a honeycomb lattice by replacing each site with a trian- 
gle, (b) Topologically equivalent representation. Sites within 
a unit cell r are labeled by a = 1, ■ • ■ , 6. The dotted line 
indicates its vertical boundary. 
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sibility (still very distant, at present) that some of this 
interesting physics might be seen on higher energy scales 
and in simpler systems than was previously realized. 

Model Hamiltonian: The Kitaev model on the triangle- 
honeycomb lattice can be written as follows: 

n = E + E + E J *rt°i 



where U*. 
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where erf (a = x, y, z) is a Pauli matrix on site i. In the 
summation, i and j are nearest neighboring sites con- 
nected by a a-link. (See Fig. [TJ) We assume that J a 
and J' a are positive in this paper (l4| . This model on a 
honeycomb lattice was studied originally by Kitaev • 
Since then, there have been some interesting develop- 
ments [ll|, [H, [H, [H, As shown by Kitaev, this 



model can be mapped to a free Majorana model with con- 
served background Z2 gauge fields in an enlarged Hilbcrt 
space. The physical states are obtained by implementing 
a set of gauge constraints. Later, Feng, Zhang and Xi- 
ang [IB} proposed a different method of solution, which 
we will adopt, based on transforming spins into Jordan- 
Wigner fermions. 

In order to implement a Jordan- Wigner transforma- 
tion [l8j], we need to order the sites. In Fig. [2(b), the 
sites can be labeled by a chain index, to, and a site in- 
dex, I along the chain. For i = (m, I) and j = (to', I'), we 
define an ordering such that i < j if to < to' or if to = ml 
and I <V . For open boundary conditions (OBC) in the x 
direction, the origin of I for this ordering is defined as the 
left edge of the system, while for periodic boundary con- 
ditions (PBC), we define the origin along an (arbitrarily 
defined) "cut," as indicated by the dotted line in Fig. 
[lib). We can define the usual Jordan- Wigner fermions: 



(Tj = (aj + dj) exp [in ^ aja.j] , 



i<j 



a] = (-l) J '(2otoj - 1), 



(2) 



(3) 



where a] is a fermion creation operator and satisfies the 



usual anticommutation relations and (— 1) J = 1 (-1) for 
black (white) sites. A notable feature of the Kitaev model 
is that, despite being two dimensional, because the only 
inter-row interactions occur on z'-links, the Hamiltonian 
expressed in terms of fermionic variables is still local. 

A further simplification occurs when we decompose 
each complex fermion into two Majarona fermions Cj and 



dj , as follows 



c = -i(a] - 
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TL can thus be rewritten as 
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fined for all pairs of sites ij connected by either a z- 
link or a z'-link. An ordering convention is adopted so 
that the products CiCj and didj are always taken with 
i < j for all links. While Eq. is exact for OBC, for 
PBC an additional boundary term must be added: TL — > 
TL + Ylcut iJ' a Fm.CiCj where the sum runs over links across 
the cut, and F m = exp[in Y.k&m a l a k] = YlkGm °k d k with 
k summed over all sites on chain m that contains link ij. 
This boundary term is crucial to obtain the correct topo- 
logical ground state degeneracy and seems not to have 
been treated in previous studies. The operators F m com- 
mute with each other and with 7i, and so are constants 
of the motion with eigenvalues F m = ± 1 . 

Manifestly, Uij is idempotent, and can be thought of 
as a background Z2 gauge connection. The set of Uij 
commute with each other, and for OBC, with the Hamil- 
tonian, as well. Thus, we can label distinct sectors of 
Hilbcrt space by the set of eigenvalues {Uij}. The Hamil- 
tonian operator in this sector is H({Uij}); it is precisely 
of the same form as in Eq. (01, but with the operators U 
replaced their eigenvalues, Uij — > Uij = ±1. This reduces 
TL({Uij}) to a quadratic form! 

For PBC which we will focus on, the situation is some- 
what more subtle since Uij on any z'-link anticommutes 
with F m for the chain m containing i or j, and con- 
sequently does not commute with Ji. In this case, we 
instead label the sectors of Hilbert space by a set of 



fluxes. 



nTj &ij where p labels the elementary pla- 



quettes (triangles and dodecagons) and the product is 
over the z and z' links surrounding the plaquette. (For 
plaquettes straddling the cut, two corresponding -F m 's 
should be included in the product, as well.) For trian- 



gle plaquettes, 



Uij. Since [H,4> p ] 



0, 



the sectors of Hilbert space can be uniquely specified by 
the eigenvalues of the local flux operators, <f> p = ±1, and 
by two global Z2 fluxes, <I> X and computed along an 
arbitrarily chosen contour encircling the system in the x 
and y directions, respectively. The Hamiltonian in each 
sector, H($ , {4> P }) , is again a quadratic form, obtained 
by replacing Uij and F m in TL by ±1 consistent with 
the specified flux pattern ($, {</> p }) with the understand- 
ing that the true eigenstates of TL should be obtained 
by applying an appropriate projection operator on the 
eigenstates of the quadratic Hamiltonian [191 ]. The cor- 
responding ground state energy is Eq($, {<fi P }). 

Broken time reversal symmetry: Under time reversal 
T, Uij — > —Uij on any z-link and z'-link; F m — > F m for 
any to. Thus <j) p is even on dodccagonal and odd on trian- 
gular plaquettes. The true ground-state is found by iden- 
tifying the sector or sectors in which Eq is minimal. In 
all cases, the ground-state of this model must break TRS, 
and must therefore be at least two- fold degenerate [20| . 
To see this, suppose a set ($, {<f> p }) which minimizes Eq. 
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There must exist a distinct time-reversed set ($ T , {0j}) 
such that E ($ T , {4%}) = £ ($,{<M) where <S T = $ 
and (f)p = <f)p ( <f>p = —(ftp) on dodecagonal (triangular) 
plaquettes. This twofold degeneracy due to TRS is in 
addition to the topological degeneracy discussed below. 

Uniform fluxes: To find the global ground-state or 
ground-states, we need to identify the sector in which 
Eq is minimal. For generic J a and J' a , we have not yet 
obtained an analytic solution to this problem. However, 
we have analyzed [2lj the problem using degenerate per- 
turbation theory, both in the limit J a <C J' a and in the 
opposite J a ^> J' a limit. Although the effective Hamil- 
tonian looks quite different in these two limits, in both 
cases it is straightforward (but tedious) to show that Eq 
is minimized by a uniform flux configuration {<fi p } = {1} 
(or the flux configuration related by T) . Numerically, we 
have computed Eq as a function of the ($, {0 P }) for a set 
of finite size systems with about 10 4 sites and for vari- 
ous set of J a and J' a . In all cases, we found that Eq is 
minimized by the same uniform flux configuration. Thus, 
the lattice translational, point group, and spin symme- 
tries are preserved. For a finite system with PBC, 
the energy difference between ground states with differ- 
ent global fluxes decays exponentially with the linear di- 
mension of the system for gapped Majorana fermions. 
In the thermodynamic limit, for the Abelian CSL, the 
ground state degeneracy is fourfold coming from four pos- 
sible values of the global fluxes, <F X = ±1 and $ y = ±1. 
However, for the non-Abelian CSL, the ground state is 
threefold degenerate because the projectio n op erator an- 

3" 



We thus 



nihilates one global flus ground state 
focus on the uniform flux (ground state) sector. The cor- 
responding free Majorana Hamiltonian, 7iijF, is of the 
same form as Eq. (j4]) with all Uij and F m replaced by 1. 

Huf can easily be exactly diagonalized. Each site is 
labeled by a unit cell index r and a site index a = 1, ■ • ■ ,6 
within the unit cell as shown in Fig. [IJb) so 

HuF = [4(c r ,4Cr+e 2 ,l) + J x(c r ,5C r ,6 + C r ,3C r ,2) 

r 

+ ^y(Cr,2Cr+ei,5) + Jy(c r .lC r .3 + C rfi C r ^) 
+ J' z {Cr^C rfi ) + J z (Cr,lCr,2 + C rA C T ^) ,(5) 

where ex,2 = x/2 ± y-\/3/2 are two unit vectors shown 
in Fig. Wlb). We define a 6-component spinor field 

-?'k-r , 



■!/> Q (k) = J2r e c r,o;/v2iV, where N is the number 
of unit cells and k lies in the first Brillouin zone. Then, 
Huf = Ek' ( /' t ( k )^( k )' ( /'(k) with H(k) a 6 x 6 hcrmitian 
matrix, and H*(k) = —H(-'k). Since c rjQ , is a Majorana 
fermion, ip a (— k) = ^(k) which means that ip satisfies 
the usual anticommutation relations: {^(k), tpp(q)} = 
<5kq<5a/3- Consequently, the problem is equivalent to a 
familiar free fermion problem, with the understanding 
that we must identify the states at k and k to avoid 
double counting of states. There are six bands, e n (k), 



corresponding to six sites per unit cell. The spectrum is 
gapped except when D(k) = Det[_ff (k)] = 0. 

Chern number and edge states: To simplify the fur- 
ther discussion, we will henceforth consider the case 
J x = J y = J z = J and J' x = J' y = J' z = J'. We note 
that these conditions are implied if we want the Hamil- 
tonian to respect the various mirror symmetries of the 
lattice pictured in Fig. [lja). 

It is straightforward to show that D(k) = only when 
k = (n, — and J' = \/2>J . Thus, there is a quan- 

tum phase transition between two fully gapped phases 
which occurs at J' = V^J, as shown in Fig. EJa). These 
two phases are topologically distinct as can be seen by 
directly computing the Chern number v: for J' > i/3J, 
v = 0, while for J' < \/3J, v = ±1. That v = ±1 for 
J' < V3J can be understood as follows: As J — > oo, the 
spectrum of Eq. ([5]) to first order in J' is gapless with 
Dirac-cones, as in graphene. Terms second order in J' / J 
generate an effective Haldane mass term which breaks 
TRS 21 1, gaps the Dirac cone, and gives v = ±1 [23j |. 
The same conclusion can be reached by computing the 
edge states of the system on a cylinder, since they re- 
flect the topological character of the bulk state [24| . In 
the small J (topologically trivial) phase, there is a gap 
about zero energy in the edge state spectrum, as well as 
the bulk. However, in the large J (topologically non- 
trivial) phase, there is an edge mode which disperses 
through zero energy. This is illustrated in Fig. [2fb) , 
in which the spectrum for J' = J of a relatively large 
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FIG. 2: (a) Bulk energy gap, A as a function of J'/J for the 
symmetric case, J a = J and J' a = J' . The quantum critical 
point is at (J'/J) c = v3. For small J, the phase is a topo- 
logically trivial CSL with Abelian excitations while for large 
J, the phase is a topologically non-trivial CSL with vortex 
excitations obeying non-Abelian statistics. The J — » oo spec- 
trum is Dirac-like, with a mass-gap that vanishes as J' 2 /J. 
There is a cusp at J'/J = 2y / 3/3, where the k with minimum 
energy jumps, (b) Energy levels calculated for J = J' on a 
cylinder with approximately 10 4 sites and with OBC along in- 
direction. Note that only half of the states are independent. 
There is one edge state (blue dotted line) crossing with zero 
energy, corresponding to Chern number v = ±1. 
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cylinder (of order 10 4 sites) is plotted as a function of k y . 
We see that there is exactly one edge state which crosses 
the zero Fermi energy, corresponding to v = ±1. Thus, 
there exists a topological quantum phase transition at 
( J'/ J) c = where the bulk gap vanishes and the edge 
states reconnect. [See Fig. HJa).] 

Non-Abelian anyons: The ground state with uniform 
{4> p } = {1} is, by definition, vortex-free. Vortex exci- 



tations are defined on plaquettes with 



= — 1. Since 



consistency requires J7 p fyp = 1 , vortices can only be cre- 
ated in pairs. In all cases that we have studied, the vortex 
excitations always cost a non- vanishing energy, i.e. they 
are gapped. It is presumably a peculiarity of the Kitacv 
model that the vortices have no dynamics, even when 
they interact with each other. 

It is to be expected that the vortex excitations in a 
state with odd Chern number obey non-Abelian braid- 
ing statistic s 1 101 , similar to the physics of p + ip super- 
conductors [25| and the Pfaffian state of 5/2-FQHE [20 |. 
Specifically, in the odd Chern number phase, each vortex 
carries a Majorana zero energy mode. A pair of Majo- 
rana zero-energy modes on two well separated vortices 
constitute a single complex fcrmionic mode with zero en- 
ergy Taking into account the four possible global fluxes 
and the effects of the projection operator [l||, the ground 
state degeneracy approaches 2 n+1 for a system of 2n well 
separated vortices. We have numerically calculated the 
spectrum for various configurations of even numbers of 
vortices in a set of finite lattices with sizes about 10 4 
sites. In all cases, we find the expected number of zero 
modes. Thus, we have confirmed that the ground state 
degeneracy for 2n vortices on a torus is 2™ +1 , independent 
of their positions or whether the vortices arc associated 
with a triangular or a dodecagonal plaquette, so long as 
the vortices are separated by a distance large compared 
to the inverse energy gap. We interpret this correct de- 
generacy as confirming the non-Abelian statistics of the 
vortices. 

In the small J (y = 0) phase, the vortices have Abelian 
braiding statistics. When J <C J', by mapping the model 
to an effective Z 2 gauge theory on a "dual" honeycomb 
lattice with spontaneously broken TRS (HI, we find the 
vortex anyons on both triangle and dodecagon plaquettes 
are bosons. Interestingly, they are relative semions. 

Discussion: As far as we know, the present model is the 
first exactly solvable model of a CSL 27[. The ground 



state spontaneously breaks TRS and is gapped, but does 
not break any spin rotational or lattice symmetries of the 
Hamiltonian. It is difficult to identify a state by what it 
does not do, which is a problem with spin liquids. Pre- 
sumably, a CSL will have signatures of TRS breaking 
such as a zero-field Kerr effect [2|| and thermal Hall con- 
ductance More generally, the effective low energy 
theory of the non-Abelian phase is probably related to 
SU (2)2 Chern-Simon theory (30| . The problem remains, 
of course, where to find this state in real materials. 
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Note added: Some of the present result were antici- 
pated in a remark made by Kitaev in Ref. [Toj . 
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